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Abstract 

General theoretical results via a Hamiltonian formulation are developed for zonal shear 
flows with the inclusion of the vortex stretching effect of the deformed free surface. These 
results include a generalization of the Flierl- Stern- Whitehead zero angular momentum the- 
orem for localized nonlinear structures (on or off the /3-plane), and sufficient conditions for 
linear and nonlinear stability in the Liapunov sense - the latter are derived via bounds on 
the equilibrium potential vorticity gradient. 
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1 Introduction 



Numerical simulations (Marcus [I]), laboratory experiments (Sommeria et al. |2], Chomaz 
et al. [3]) and analytical treatments (Shivamoggi [4]) indicated that, in a quasi-geostrophic 
flow, the existence of a strong zonal-flow shear can lead to the formation of a strong persistent 
coherent vortex. On the other hand, large scale wave motions, or meanders, occurring in the 
oceanic circulations are believed to tbe caused by the instability of the zonal shear flow (Kuo 
[5]). On this premise, the linear stability of a barotropic zonal shear flow render the influence 
of differential rotation (as represented by the Rossby parameter (3) was addressed by Kuo 
[6], Lipps [7j, Pedlosky [8] and [9], Howard and Drazin [ID] , Drazin et al. [11], Hickernell 
Gnevyshev and Shrira |13j . and others. These investigations were generalized (Ripa 
Shivamoggi and Rollins [J5]) to include the vortex-stretching effect of the deformed free 
surface. In this paper, we give a Hamiltonian formulation for this problem, and use this 
framework to give - 

1. a generalization of the Flierl-Stern- Whitehead [16] zero angular momentum theorem for 
localized nonlinear structures (on or off the /3-plane); 

2. a sufficient condition for linear stability in the Liapunov sense for steady states of this 
system that complements the previous result of Ripa 



3. sufficient conditions for nonlinear stability in the Liapunov sense via bounds on the 
equilibrium potential vorticity gradient. 

2 Hamiltonian Formulation 

Consider a two-dimensional quasi-geostrophic flow in which the baroclinic effects are pro- 
duced by the deformed free surface of the ocean. The governing equation (in appropriate 
units) is (Charney [17] 

§ + [<M = o (l) 

where <f> is the stream function and q is the potential vorticity - 

^V 2 0-0 + / (2) 



f being the Coriolis parameter and, 



[4 a] se— — -— — (3) 

dx dy dy dx ' 



The Hamiltonian for this system is 



H = - l - / fad A. (4) 



A 



assuming that (p = on the boundary C enclosing the area A where the flow occurs. 
Take q to be the canonical variable and the skew-symmetric operator J to be 

J = -[«,(•)]■ (5) 
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The Hamilton equation is then given by 



which is just equation (1), as required! 

The Casimir invariants for this system are the solutions of 



SC 

which implies 

SC 



SC 



(7) 



F'(q) (8) 

uq 

or 



C = jF(q)dA. (9) 



^ = ^ do) 



or 



C may be understood by noting that q is actually a Lagrange invariant, so any integral 
function of q would be conserved! 

3 The Beltrami States 

The Beltrami state is the minimizer of H keeping C constant, and is given by 

5H _ X 5C 
Sq Sq 

-4> = \F\q) (11) 

or 

q = P{<f>). (12) 

4 Trans lational Symmetry Solutions 

If this system has translational symmetry along the x-direction, the latter is characterized 
by an invariant M given by 

which implies 

M= yqdA. (14) 

J A 

On the other hand, existence of a traveling wave solution for this system implies that 
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q(x,y,t) =q(£,y),£ = x-ut (15) 

which means, in turn, 

dq dq 

- U d~x = W (16) 
On using equations (6) and (13), equation (16) implies 

Uj lq~ = ^ (17) 



or 



J-(H-uM) 
dq 



q^{H-uM) 



On using (4) and (14), equation (18) leads to 

[q,<P-uy]=0. (19) 

5 Generalization of Flierl-Stern- Whitehead Zero Angular Momen- 
tum Theorem 

Multiply equation (19) through by £, and integrate over A; this leads to 

/ [q,(j ) -uy^dA = 0. (20) 

J A 

Assuming localized traveling structures, (20) reduces to 

[((f)- uy)q y dA = (21) 

J A 

or 

/ cj)fdA + u I qdA = 0. (22) 

J A J A 

Assume Beltrami states given by 

q = acf) (23) 

a being a constant, which is a special case of (12) and is of added interest because it then 
yields an H which is quadratic in cp\ The Beltrami condition is equivalent to the "slow 
variation" restriction used by Flierl, Stern, and Whitehead (FSW) [16] and quantifies it 
precisely 

On using (23), equation (22) becomes 

/ {au + f)4>dA = (24) 

J A 



4 



which is a generalization of the FSW zero angular momentum theorem for localized non- 
linear structures that results in the /3-plane approximation (/' = const = j3)- 



(fxLA = 0. 



(25) 



Note that in the generic case, more general localized nonlinear structures with non-zero 
angular momentum are possible, especially off the /3-plane! On the other hand, even on the 
/3-plane, observe that the traveling localized nonlinear structures are not restricted by the 
FSW theorem if they satisfy the critical condition - 



u = . 

a 



(26) 



This condition, it is to be noted, is not restrictive because the constant a is arbitrary. It 
may be noted that the present formulation could be trivially generalized to include zonal 
shear flows. 

6 Linear Stability in the Liapunov Sense 

The stationary state of the system in question is a stationary point of the functional - 



I[q]=H[q]+([q] 



+ F(q) 



dA 



and is given by the vanishing of the first variation of I[q] 



8I[q) 



+ F'(q)]6q dA = 



from which, 



cj) = F'(q) orq = P(<f>). 
Next, the second variation of I[q] is given by 



6 2 I[q] 



i(V50) 2 + \{54>f + 1 -F"(q)(5q) 2 ] 



dA. 



(27) 



(28) 



(29) 



(30) 



A sufficient condition for linear stability in the Liapunov sense (Arnol'd [18J) is that 8 2 I 
is positive (or negative) definite for all suitably smooth variations 5q. This requires, from 
(29), that 



F"{q) > 0, V(x,y) G A 



(31) 



!The net angular momentum is 

L = JJ A (xv — yu)dxdy 

= - J f A (x(j) x + y4> y )dxdy 

= ff A 2(pdxdy. 
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which, on using equation (28), implies 



T q * < 32 > 

or 

(l)(f)*°- v( -» e - « 33> 

If the system has translational symmetry along the x-direction, then (32) becomes 

dq\ ( d .— 



For the parallel zonal shear flows, u = u(y) , the linear stability criterion (33) becomes 

(u" - (3 -u)(u- c) > (35) 

with the equality (where the potential vorticity attains its maximum) valid at y = y c where 
u — c- in agreement with the normal-mode result (Shivamoggi and Rollins [15]). In the 
absence of the planetary rotation, this result reduces to Fjortoft's theorem. 

The above result was also given by Ripa [H] by using an energy method that made 
use of invariants associated with explicit symmetries of the system. It may be noted that 
the present development complements this formulation by using instead Casimir invariants 
associated with the non-canonical nature of the Eulerian formulation of the system. 

Note that the stability criterion (34) is applicable to an arbitrary perturbation q which 
allows I[q] to exist. Thus the normal- mode result is generalizable to arbitrary perturbations 
and implies even the normed stability: 

THEOREM: The steady state (0, q) is linearly stable in the Liapunov sense with respect 
to the perturbation norm - 

\\(5<f>,6q)\\ = {6 2 I[q]y (36) 
if the Casimir function F(q) satisfies 

F"(q)>0, V(x,y) G A. (31) 

7 Nonlinear Stability 

Nonlinear stability theory seeks to establish stability of a steady state by shown that it is a 
local extremum of energy with respect to /imte-amplitude, Casimir-preserving perturbations. 
In order to address nonlinear stability, one needs to examine the sign of the exact invariant 
(Arnol'd [18J) - 

Al[q] = I[q] - I[q] (37) 

instead of \o~ 2 I which was considered in the linear stability problem. If AI is positive 
(or negative) definite, then it can be used to define a perturbation norm, and the steady 
state (29) is nonlmearly stable in this norm. This development needs additional convexity 



G 



conditions to be imposed on the functional I[q] in order to insure that the latter is convex 
in a small but finite neighborhood of the steady state (29). 
Let us write AI[q] as follows - 

AI[q] = I[q + q]-I[q] (38) 

where q is the /mite-amplitude perturbation imposed on the stationary state q. 
Using (27) and (29), we have, 

AI[q] = - J \UdA + J [F(q + q)- F(q) - F\q)q]dA. (39) 

Note, 

AI[q] = 5 2 I + O(0 3 , q 3 , 4> 2 q, 4>q 2 ) (40) 

so, 5 2 I is only the sma//-amplitude approximation to All 
Suppose F(q) satisfies the convexity condition - 

< ai < F"{q) < /3i < oo, Vg (41) 
where cti and f3\ are real numbers. So, 

l - ai q 2 < [F(q + q)- F(q) - F'(q)q] < \^q 2 (42) 
and hence, from (39), we have 

I J[^ q + ai f ]dA < A /[g] < I jf [-^ + /? 1( f]dA (43) 

which shows that the convexity condition (41) is sufficient to establish the positive (or nega- 
tive) definiteness of the exact invariant AI[q]. Further, (43) implies that a perturbation norm 
defined by 



\q\\ = {AI*[q}}2 (44) 



where, 



Ar[ q ] = ±Jj-tq + ai q 2 ]dA (45) 

is bounded from above. 

THEOREM: Suppose that the Casimir function F(q) satisfies the following convexity 
condition - 

< ai < F"{q) < fa < oo, Vg (41) 
for some real constants a.\ and (3\. Then, the stationary state q, determined by 

= F'(q) (29) 
is nonlinearly stable in the Liapunov sense with respect to the perturbation norm: 
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\\q\\ = {AI*[q}}^. (44) 

(41) implies, it may be noted, explicit bounds on the equilibrium potential vorticity 
gradient which signify conditions stronger than (31) that is required for linear stability. For 
nonlinear stability, the potential vorticity gradient not only needs to be positive (as for linear 
stability) but it must also be bounded from above. 

8 Discussion 

In this paper, a Hamiltonian formulation has been used to develop general theoretical results 
for zonal shear flows with the inclusion of the vortex stretching effect of the deformed free 
surface. One such result is a generalization of the FSW zero angular momentum theorem 
for localized nonlinear structures (on or off the (3- plane). Off the /3-plane, more general 
localized nonlinear structures with non-zero angular momentum are shown to be possible. 
On the other hand, even on the /3-plane, localized nonlinear structures are shown not to be 
restricted by the FSW theorem provided they propagate at certain critical speeds. Another 
general theoretical result is concerned with sufficient conditions for linear and nonlinear 
stability in the Liapunov sense - the latter have been derived via explicit convexity conditions 
for the Casimir function and hence bounds on the equilibrium potential vorticity gradient. 
For nonlinear stability, the latter not only needs to be positive (as for linear stability) but 
it must also be bounded from above. 
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